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Abstract 

Based on the finite-temperature AdS/CFT correspondence, we calculate the ratio of 
shear viscosity to entropy density in any Lovelock theories to any order. Our result 
shows that any Lovelock correction terms except the Gauss-Bonnet term have no 
contribution to the value of rj/s. This result is consistent with that of Brustein and 
Medved's prediction. 

PACS: number(s): 11.25.Tq; 04.50.-h; 04.70. Dy; 11.25.Hf 

Keywords: AdS/CFT, KSS bound, Lovelock theory, AdS black brane. 



Stimulated by the conjecture of AdS/CFT correspondence [IH3], string theory has attracted 
a lot of attention, especially after the discovery that some theoretical results of the dual 
theory are consistent with that of the RHIC experiment, say, the ratio of the viscosity to the 
entropy density [HE]. Recently, it was conjectured, based on the AdS/CFT correspondence, 
that for all possible nonrelativistic fluids, there may exist a universal lower bound (the KSS 
bound) on the viscosity/entropy-density ratio (we set G = c = h = k B = 1) [B] 

s 4i 

This bound received great supports from several kinds of field theories [IHS], as well as 
the case with chemical potential in the theory |10[ITT]. However, more recent work on the 
higher derivative gravity theories (see [i~2rli~8] ) showed that the KSS bound is violated 
when the dual gravity is enlarged to include a stringy correction (see [19] for more about 
the KSS bound in higher derivative gravity). This correction is frequently referred to as 
the quantum correction, since in CFT side this is a correction of the 't Hooft coupling 
A = g\ M N c . It is of particular significance to consider the 1/A correction when we are 
dealing with non-extremely strong coupling fluids. 
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Recently, the authors of [20] predicted that all Lovelock terms higher than the second 
order(the Gauss-Bonnet term) do NOT contribute to the value of 77/s at all, and this 
prediction was partially confirmed in [21 J for the third-order Lovelock gravity. In this 
paper we calculate the viscosity/entropy-density ratio directly in the Lovelock theory to 
any order, trying to make a complete verification of the prediction, and indeed, our result 
provides a direct support of this prediction as will see below. 

We start with the Lovelock theory of gravity. This is one of the most general second order 
gravity theories in higher dimensional spacetimes and is free of ghost when expanding on 
a flat space [22J and hence is of particular interest. The Lagrangian density for general 
Lovelock gravity in D dimensions is C = Sm=o c m£ m , where C m is given by [2"3] 

&m C* Pl^l D Pm^m ( r )\ 

c m is the m'th order coupling constant, [D/2] denotes the integer value of D/2 and the 
Greek indices A, p, a and k go from to D — 1. The symbol R\ a pK is the Riemann tensor 
in D-dimensions and ^Ik{..-p™kZ ^ s the generalized totally antisymmetric Kronecker delta. 
The term £ — V~ 9 * s the cosmological term, while C\ = y/— fi"^ 1 ^ 1 R\ iai PlKl /2 is the 
Einstein term. In general C m is the Euler class of a 2m dimensional manifold. 

Variation of the Lagrangian with respect to the metric yields the Lovelock equation of 
motion 
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As is shown in [21] , there exist static exact solutions of Lovelock equation. Let us consider 
the following metric 



ds 2 = -f(r)dt 2 + —— + r 2 Yl lijd^dx\ (4) 

/ \ ) i,j 

where r )ijdx l dx^ represents the line element of a (D — 2)-dimensional Einstein space. With 
this ansatz, we have 

Kijki = ^iliklji ~ luljk), Kij = k{D - 3)7y, TZ = k(D - 2)(D - 3). (5) 

where k is the curvature constant, whose value determines the geometry of the horizon. 
Without loss of the generality, one may take k — 1, 0, or — 1 representing sphere, flat and 
hyperbolic respectively. 

Using this metric ansatz, we can calculate Riemann tensor components as 
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Rtr tr = - y> R * j = R ™ 3 = ~Y r 5 ^ R v kl = i 6 ^ 1 - 5i ^ h ) ■ (6) 

Substituting (J5J) into (J3J) derives a simple equation 



W[4>\ = £ *»V> m = r5~i' (7) 

m=0 ' 

where ■?/> = r _2 (/t — /), /i > is a constant of integration which is related to the ADM mass 
by 

where Vd-2 is the volume of the (D — 2)-dimensional hypersurface and Gd is the Newton 
constant. In (J7J), we also defined c m = ^^2m-iy. Cm an< ^ n ^ s an m ^ e g er with < n < [D/2\. 
In this paper we are considering AdS black brane in Lovelock gravity, so we have cq = — 2A 
with the cosmological constant A = —(D — 1)(D — 2)/2/ 2 and c\ = 1. 

We would like to extract some information from the Lovelock black brane, such as their 
thermodynamic properties. One quantity which is of particular interest is the entropy S. 
Generally speaking, one can obtain the entropy of a black hole in higher derivative theories 
by using the thermodynamic relation S = —dF/dT with F the free energy and T the 
Hawking temperature. By doing so one finds that the entropy of the Lovelock black brane 
is given by [25] 

_ A m(D-2) 

where r + is the event horizon of the black brane which is the positive root of f(r + ) = 0. In 
the present paper, we mainly focus on the case where k — 0. In this case we have a simple 
formula for the entropy density of the Lovelock black brane 



r D-2 

1+ 
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(10) 



and now, r + is a solution of 4>{ r +) = 0- The Hawking temperature of this case is given by 

T-^r + . (U) 

In what following, we would like to see the waves generated by a metric perturbation of the 
background. Generally speaking, there are scalar, vector and tensor modes depending on 
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the rotation symmetry. In this paper, we only study tensor perturbations which is closely 
related to the shear viscosity as will see below. 

We now add a small tensor perturbations to the solution (j4j) 



Sgij = r 2 (fi(t,r)hij(x' 1 ) , (others) = 



'121 



where 4>(t, r) represents the dynamical degrees of freedom. Here, hij are defined by 



V k V k hij = k 2 h 



V l h {j = , 



l ij h i3 = 0. 



(13) 



Here, V 1 denotes a covariant derivative with respect to 7^ and k 2 is the eigenvalue playing 
a role of momentum. 

With these definition, one can obtain the first order perturbation equation of the Lovelock 
equation ([3]) [26] 
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where 5R a b cd represents the first order variation of the Riemann tensor and we have intro- 
duced a new quantity a m = mc m {m > 0). As shown in [26], it is straightforward once we 
know the expressions of quantities SRtfi , 5R r i'^ and 5Rij kl . Then the calculation becomes 
a mathematical game and the result is ready-made [26] 
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Using the Fourier decomposition 



(17) 



we obtain the linearized equation of motion for (f>(r): 
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It is convenient to introduce a new dimensionless coordinate u = (r + /r)( D ~ 1 < ) / 2 with r + 
the event horizon of the black brane. In this coordinate frame, u = corresponds to the 
boundary and u = 1 the horizon. The linearized equation of motion f ll8p then becomes (for 

K = 0) 

, q'(u) i// \ u) 2 / \ D — l h'k 2 ,. . . . 

0" (u + (u) + 2g _ 6 - > m - • p _ 8 - #u = (19) 
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where 



^(u) = -rl~ u ilj(u)h(u)u—, (20) 
2 2 

— i~\ — U) fc =m — r\ — ^' ( 21 ) 
[D — l)r + (L> — l)r+ 

and the prime denotes the derivative with respect to u. 

Now we shall calculate the shear viscosity in Lovelock gravity theories. Generally speaking, 
the shear viscosity 77 can be calculated via Kubo formula, 

, Im(G H (w,0)) , . 

77 = - lim — i — Ljl , 22 

where G R is the retarded Green's function 

G R (u,k) = -i I dtdxe- i% - s 6{t) < [6(x)6(0)} >, (23) 



with O some boundary CFT operators. According to AdS / CFT correspondence, the Green's 
function can be calculated from the dual gravity side via the Gubser-Klebanov-Polyakov/Witten 
relation [2J|3] 

where <fi is the bulk field and 0o is its value at the boundary, i.e., 0o = lvoa u ->o 4>{ u ) • 
Extracting the part of S c i that is quadratic in and inserting the solution of the linearized 
field equation we may get a surface term in four dimensions by using the equation of 
motion, 



Sd[M =J ^ £> _ 1 o (-fc)<2(fc,^)0o(fc) 



(24) 



M=0 

where u = (r + /V) as defined previously. In this way, we obtain the following relation 

for the retarded Green's function [27] 

(25) 

u=0 



G R (k) = 2G(k jU ) 
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where the incoming boundary condition at the horizon is imposed. The shear viscosity then 
can be calculated by using (|22|) . 



In the following we would like to calculate the shear viscosity, following the procedures 
introduced above. The main task is to solve the equation of motion f)19p in hydrodynamic 
regime i.e., small uj and k. To solve the wave equation (TT9|) we first examine the behavior 
around the horizon where u — 1. For this purpose it is convenient to impose a solution as 

^u) = (l-uyF(u), (26) 

with F{u) regular at the horizon. Substituting (126]) into the wave equation f lT9|) and leaving 
the most divergent terms, we can obtain 

» = (27) 

where we have used the relations 



g(u -> 1) = -g'(l)(l -u) + 0((1 - u) 2 ), (28) 
ip(u 1) = -^'(1)(1 -u) + 0((1 - uf). (29) 

In present paper we choose "— " sign in eq. f )27|) for convenience. 

To get the viscosity via Kubo formula fl22|) . the standard procedure is to consider series 
expansion of the solution in terms of frequencies up to the linear order of u, 

F{u) = F (u) + vFi(u) + 0(v 2 , k 2 ). (30) 

Then the equation of motion f[T9l becomes the following form up to 0(u), 

[g(u)F'(u))' - v (j^9(u))' F(u) - ^-g( u )F\u) = 0. (31) 

After substituting the series expansion ( 1301) into the equation (13T1) . we obtain the following 
equations of motion for F Q (u) and F\{u) 



[g(u)F^u)}' = 0, (32) 
[g(u)F[(u)}' - (-^(«))' F o(u) = 0. (33) 

By requiring that the functions Fq(u) and Fi(u) are regular at the horizon one gets the 
following results 



F[{u) 




(34) 
(35) 
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where again we have used the relation ([28]) and the constant C can be determined in terms 
of boundary value of the field, i.e., C — 0o( 1 + ^{v) ) . 



Now we shall calculate the retarded Green's function. Using the equation of motion, the 
action reduces to the surface terms. The relevant part is given by 
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Near the boundary u = e, using the perturbative solution of <f>(u), we get 



(36) 



n=0 



( j ) \e) = v g -^±^ + 0{v\k 2 ) 

-m^ o+0i ^ k2) - (37) 

Therefore we can read off the correlation function from the relation (125*1) , 

G >- fc) =-iftk(ST) J ' (1)+0(w2 '* 2) ' (38) 

where contact terms are subtracted. Then the shear viscosity can be obtained by using 
Kubo formula (J22D, 

" - -isk (m) g,{1) - (39) 

The ratio of the shear viscosity to the entropy density is concluded as 

i = (40) 

s 4vr^'(l) 

From f[2Ul) we have a relation g'(l) = — r^r D ip'(l)h(l) and h(l) can be obtained from ffTB"]) 
by inserting k = 

/i(l)=r^ 4 (l- (D-l) (£-4)5(^2 
It is straightforward to show that 

where we have defined A = c 2 . This result is exactly the one predicted in [2D] . 

In summary, we have computed the ratio of shear viscosity to entropy density for any 
Lovelock theories. Our result shows that any correction terms except the Gauss-Bonnet 
term do not affect the value of rj/s, and this confirms the prediction made by [20]. During 
our calculation, we have chosen a vanishing curvature constant k. Actually, our result is 
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still valid (for leading term) for nonzero k if we focus on a large black brane. In the large 
black brane limit, both the entropy and the viscosity have the same leading terms as those 
of k = 0. This can be seen by noting the expressions of entropy density and viscosity. From 
([2]), the entropy density of the Lovelock black brane with nonzero k can be expanded, in 
the large black brane limit (i.e., <C 1), to the first order as 



With the same spirit one can also expand the shear viscosity to the first order in the large 
black brane limit. This can be done by repeating the previous procedures and noting that 
h(l) = h(u = 1) can be obtained from (1161) . In this way, the shear viscosity for nonvanishing 
k can be expanded to the first order as 



From ( l4"2l) and ( 143]) it is obvious that the leading terms of the entropy density and the 
viscosity for k ^ are the same as those of k = 0. In other words, in the large black 
brane limit, the curvature constant k has no contribution to the shear viscosity to entropy 
density ratio for the leading term. The sub-leading terms, however, receive contributions 
from k. 

So far we are confident with the violation of the KSS bound while we are not sure the 
existence of a universal lower bound of r]/s. A great progress alone this line appeared 
several months ago when the authors of (28] gave a proof for the existence of a universal 
bound of rj/s for any ghost-free extension of Einstein theory. However, the work in [T4] 
shows that the causality violation of the dual gauge theory may put constraints on the 
coefficients of higher derivative terms and this in turn will put constraints on the value 
of rj/s. Then it is natural to ask if the lower bound still exists as these constraints are 
taken into account. Recent progress made by Camanho and Edelstein in [29] provides us 
with an answer that the causality violation, as expected, may impose a constraint on the 
bound of the r)/s at least for cubic Lovelock gravity. For completeness, we briefly catch 
some important results from [29], so as to compare the lower limit on rj/s from causality 
violation and the result in the present paper. Actually, the formula of the ratio between 
rj and s obtained in [29J is not different from our result (j4"T|) . What is new of their result 
is that by imposing a condition so that the causality violation can be avoided, they found 
constraints on the coefficient A (or C2 as defined) in (T4T1) . For any order (n > 2) Lovelock 
gravity, the condition to be free of causality is that 



+ 




(42) 




(43) 




(44) 



S 



where 7 = —2, —1,2/(D — 4) represent helicity zero, helicity one and helicity two gravi- 
ton, respectively. Therefore, though any correction terms higher than the second order of 
Lovelock gravity do not manifestly contribute to the ratio of viscosity to entropy density, 
it does not mean that they are irrelevant to this ratio. Through ( |44p we see these terms 
impose a constraint on the value of C2 (or A) thus in turn affecting the lower bound for 
r)/s. 
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